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Abstract 

We apply self consistent field theory to twist 
grain boundaries of block copolymer melts. The 
distribution of monomers throughout the grain 
boundary is obtained as well as the grain bound- 
ary free energy per unit area as a function of 
twist angle. We define an intermaterial dividing 
surface in order to compare it with minimal sur- 
faces which have been proposed. Our calculation 
shows that the dividing surface is not a minimal 
one, but the linear stack of dislocations seems to 
be a better representation of it for most angles 
than is Scherck's first surface. 



1 Introduction 

Bulk equilibrium properties of diblock copoly- 
mer melts are relatively well understood [jjj. In- 
compatibility of the monomers comprising the 
two blocks drives the system toward ordered 
structures in which the number of contacts be- 
tween dissimilar monomers is reduced, subject 
to various constraints. These ordered phases, 



of which the simplest is lamellar, are thermody- 
namically stable below some order-disorder tran- 
sition temperature. 

When the system is taken below this tempera- 
ture, the lamellar phase is nucleated typically in 
distinct grains which differ, one from the other, 
by the orientation of the lamellae within them. 
The interface between lamellae of different grains 
constitutes a grain boundary, which can be con- 
sidered an equilibrium structure arising from a 
constraint that imposes the different orientations 
of the lamellae of the two grains. 

Because the lamellae of block copolymers lack 
any internal order, their grain boundaries are 
simpler than those between grains of crystalline 
solids. While the latter are combinations of five 
different independent boundaries, the former can 
be decomposed into only two independent ones. 
In the kink grain boundary, the normals perpen- 
dicular to the lamellae of the two grains define 
a plane which is perpendicular to the plane of 
the boundary. Kink grain boundaries have been 
studied recently both experimentally [||, ||, |j and 
theoretically || §, @j. In the twist grain bound- 
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Figure 1: Schematic of a twist grain boundary 
in which one looks down the x axis. As x — > oo, 
the IMDS approach that of the lamellar bulk, 
shown in dashed lines, while as x — * — oo, the 
IMDS approach that of the lamellar bulk shown 
in solid lines. Normals to the planes of the bulk 
lamellae make angles of ±a/2 with respect to 
the y axis. 



ary (TGB), the plane defined by the normals is 
parallel to the plane of the boundary. The angle, 
a, between the normals defines the twist angle. 
The geometry is shown in Fig. 1 where we fix the 
convention we use throughout: x is the direction 
perpendicular to the grain boundary and y and z 
are in the plane of it; normals to the lamellae of 
the two grains are at angles ±a/2 with respect 
to the y axis. Twist grain boundaries have been 
the subject of a thorough experimental study j^, 
but have received less theoretical attention. 

Two different kinds of twist grain boundaries 
have been observed. The simplest consists of a 
stack of planes. In each plane, the orientation 
of lamellae differs slightly from those in planes 
above and below it. The structure is periodic 
along the stacking direction. It is observed for 
small twist angles (a < 15°) only. A heuristic 



model for this boundary was proposed by Gido 
et al. H and checked against their experimental 
data, with good results. 

The other structure, observed for all angles, 
is quite different: it is doubly periodic. A sur- 
face which displays this double periodicity is 
Scherck's first surface, one which is characterized 
by zero mean curvature everywhere. It is shown 
in Fig. 2 for a twist angle a = 0.5 radians. This 
minimal surface was first proposed as a model 
of the TGB by Thomas et a/.|J. The reasoning 
is as follows. Suppose that one can ignore the 
components of the system, the block copolymers 
and the constituent monomers, A and B of the 
two blocks, and focus instead on the internal in- 
terfaces which divide the A-rich lamellae from 
the i3-rich ones. Suppose further that one can 
ignore the finite thickness of these internal in- 
terfaces and approximate them by a suitably de- 
fined surface, the intermaterial dividing surface, 
(IMDS). The area of this surface is extensive, i.e. 
proportional to the volume of the system. Un- 
der these assumptions, the bulk free energy of 
the system should be expressable as the energy 
of this surface. Minimization of this free energy 
leads to the surface with minimum area subject 
to the constraint that it separate regions of cer- 
tain volume. This is a surface of constant mean 
curvature. In particular for a symmetric system 
in which the volumes are equal, and that is the 
case for a diblock with equal volumes of A and B 
monomers, the constant value of the mean cur- 
vature is zero. Such surfaces are called minimal 
surfaces. The choice of the appropriate minimal 
surface depends upon boundary conditions and 



expected symmetries. Scherck's first surface[lC] 
was chosen in Ref || for comparison with ex- 
perimental results because it connects two sets 
of parallel planes, the normals of which form an 
angle a. The explicit expression for the surface 
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is quite simple: 
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The equation denning the surface pertains only 
in the regions for which the right hand side is 
positive. This condition defines a chessboard- 
like domain in the y-z plane (as a consequence 
of which the surface is sometimes referred to as 
"Scherck's doubly periodic surface"). Far from 
the grain boundary, whose center is at x = 0, 
the nature of the surface is easily inferred from 
Eq. |l|: when x — > — oo the r.h.s. must vanish, 
and this defines a set of parallel planes with unit 
normal (0, cos (a/2), sin(a/2)) and spacing D/2. 
The period of the IMDS is half that of the lamel- 
lar phase, D, hence our convention. Similarly, 
when x — > oo the r.h.s. must diverge, and this 
defines a set of parallel planes with unit normal 
(0, cos(a/2), — sin(a/2)) and spacing D/2. 

A different description of the grain boundary, 
due to Renn and Lubensky [|ll|] , is that of a linear 
stack of dislocations (LSD). The model of a sin- 
gle dislocation is again a minimal surface, the he- 
licoid. The whole structure arises from the stack- 
ing of infinitely many dislocations along a line. 
In this case the dislocations have vorticity along 
the y axis with a pitch D, and are stacked along 
the z axis with a separation D/(2sin(a/2)). It 
was later shown |12| that this approach is ac- 
tually equivalent to the description employing 
Scherck's surface, up to a cos (a/2) dilation of 
the x axis, the LSD being more "compressed" 
than Scherck's surface. 

The above approaches are valuable in provid- 
ing simple models of the grain boundary which 
can be compared with experiment. They suf- 
fer, however, from the approximations which are 




Figure 2: Scherck's first surface is shown for a 
twist angle of a = 0.5 and a scale D = 1. 



inherent in the approach. Most importantly in 
the system of block copolymers, they ignore the 
physical constraints of incompressibility which 
causes the chains to stretch in order to fill the 
available volume. As a consequence, the IMDS is 
not a surface of constant mean curvature, a point 



made compellingly by Matsen and Bates [ 13 1, and 



confirmed in experiment [ 14 1. 

Therefore we have studied the twist grain 
boundary in block copolymer systems using the 
self-consistent field theory in Fourier space which 
has been successful in describing lyotropic phases 
of block copolymer melts [15]. We follow the ap- 
proach of Matsen [|j who adapted it to kink grain 
boundaries. First we will introduce the partic- 
ular implementation of the theory which is use- 
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Figure 3: Schematic of the system, periodic in 
all three directions, on which we perform our 
calculations. The original system of interest is 
recovered on letting A increase without limit. 

ful in this case. We then present results for the 
TGB obtained within this framework and com- 
pare them with Scherck's first surface and the 
LSD. We conclude with a brief discussion. 

2 Theory 

We consider an incompressible melt of n AB di- 
block copolymers each composed of N segments 
of volume 1/poj the volume of the system is 
V = nN/pQ. The polymers are modeled as Gaus- 
sian random walks with statistical step length a. 
The natural length scale of the problem is the 
end to end mean distance, aN 1 ^ 2 . To describe 
the incompatibility between A and B monomers 
the standard Flory-Huggins parameter, x> is in- 
troduced; the product xiV sets the units of en- 
ergy and temperature. 

We utilize the SCFT method expressed in 



Fourier space as in Ref. [15], a method suited 
to the study of periodic phases. A system con- 
taining a twist grain boundary, however, is only 
periodic in the coordinates y and z, but not in x. 
In order to circumvent this difficulty, we adopt 
the same strategy employed by Matsen in his 
study of the kink grain boundary Q: to express 
the desired system as one periodic in all three di- 
rections in the limit in which one period becomes 
infinite. We therefore consider the system shown 
in Fig. 3 which, in addition to being periodic in 
y and z is periodic in x with period A. We im- 
pose a reflection symmetry around x = A/2, i.e., 
(x, y) —* (A — x, y). The desired grain boundary 
free energy per unit area, T, of the original sys- 
tem is obtained from the free energy, F(A), of 
our system according to 



T = lim [F(A) - F b ] 



V 



(2) 



where the bulk free energy Fy, = limA^oo F(A), 
and the area of the grain boundary is V/A. 
As the natural unit of area per polymer is 
V/naN 1 / 2 , the natural unit of the grain bound- 
ary free energy is k-gTnaN 1 / 2 /V . Thus the 
boundary free energy can be written 



V 



kBTnaN 1 / 2 
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A^oo 



F(A) 



A 



k B Tn aN 1 / 2 



(3) 

As the system is now periodic in all direc- 
tions, we can expand all functions of position 
into a complete, orthonormal, set of eigenfunc- 
tions, fimn, of the Laplacian operator pq], eigen- 
functions which explicitly express the symme- 
tries of the system. It is clear that it is invariant 
under a rotation of tt about the x axis, i. e. under 
(x,y,z) — > (x, — y,— z). It is also invariant un- 
der rotations of it about the y and z axes as well. 
Therefore the system is invariant with respect to 
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the change in sign of any two of the coordinates. 
These considerations, together with the imposed 
reflection symmetry around x = A/2, lead to the 
choice of functions 

flmnipi ViZ) = 

cic m c n cos(lk x x) cos(mk y y) cos(nk z z), 

if I is even 
QC m c n sixi(lk x x) sm(mk y y) sin(nk z z), 

if I is odd, 

(4) 

where k x = vr/A, k y = 2ircos(a/2)/D and 
k z = 2ir sin(a/2)/D. Again D is the bulk lamel- 
lar period. 

It should be clear from Fig. 1 that the natural 
coordinates in which to express the periodicity of 
the system are (y/D) cos(a/2) ± (z/D) sin(a/2). 
Indeed Scherck's surface, Eq. |], is expressed in 
them. Translating an expansion in those coordi- 
nates into an expansion in x, y, and z, one sees 
that the parity of m and n above must be the 
same. Finally the cj are determined by orthonor- 
mality: 

V J fl mn m^)dv = 5 lv 8 mm ,5 nn , (5) 

Thus Co = 1 and q = v2 for i > 0. This com- 
pletes the specification of the basis functions. 

One might think it necessary to impose an 
additional invariance: that the calculated free 
energy of the grain boundary with twist angle 
a < tt/2 be identical to that with angle tt — a 
because, after an interchange of y and z, the 
one boundary is identical to the other. However 
the grain boundary free energy we calculate al- 
ready displays the symmetry T(a) = T(ir — a) 
without further restriction of the basis func- 
tions. This can be seen from the fact that under 
a — > tt — a, the wavevector components k y and 
k z interchange. Thus an interchange of the coor- 
dinates y and z and a relabelling of the dummy 



indices m and n suffices to make the expressions 
for the free energies of the two grain boundaries 
identical. 

We expand all functions of position in terms of 
the above basis functions. Of course the infinite 
set must be truncated in a numerical calcula- 
tion, and our computer resources impose a max- 
imum slightly below 400 functions. The results 
presented below are obtained for the choice of 5 
values for m, 5 for n and 15 for I as the results 
are more sensitive to the number of components 
in the x direction than in y or z. To be consis- 
tent with this choice, we take the corresponding 
bulk lamellar phase to be that obtained from 5 
Fourier components. 

We have chosen \N = 15, an intermediate 
value for which an intermaterial dividing surface 
is well delineated, but not so sharp as to require 
a large number of basis functions to describe. A 
value of A ~ 5aA^ 1 / 2 is sufficient typically for 
the free energy to become insensitive to further 
increases in this parameter. With our 375 ba- 
sis functions, our results for the free energy at 
a = tt/2 are accurate to within 1%. Larger val- 
ues of X-ZV would require additional basis func- 
tions. Smaller values of %iV, nearer the order 
disorder transition temperature of %iV = 10.49, 
cause periodic modulations of the dividing sur- 
face to appear which extend away from the grain 
boundary. This behavior, similar to that re- 
ported for the kink grain boundary || [|, is 
likely to be strongly modified by fluctuation ef- 
fects which are absent in the SOFT. 

3 Results 

In Fig. 4 we show results for the twist angle 
a = 0.4 w 22.9°. We have plotted, for several 
values of x, contours of constant order parame- 
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Figure 4: Slices through the grain boundary 
showing the monomer distribution at different 
distances: (a) In the limit of x — > oo, i.e. the 
bulk, (b) x = O^aN 1 / 2 , (c) x = 0, i.e. at the 
grain boundary itself. In these gray scale plots, 
the maximum absolute value of the order param- 
eter is 0.88. 

ter, the difference between the volume fractions 
of the two monomers. In these gray scale plots, 
the maximum absolute value of the order param- 
eter is 0.88. Fig. 4(a) shows a slice at infinitely 
large x, that is, a cross section through the bulk 
system. Figs. 4(b) and (c) show slices at the 
values of x = 0.2 aN 1 / 2 , and 0, that is, at the 
grain boundary itself. 

We would like to compare the distribution 
of monomers obtained in our solution with 
Scherck's first surface, which is a model for the 
intermaterial dividing surface. One way to do 
this is to calculate within our solution the value 
of the order parameter, 5(j)(r) = 4>a{y) — 1/2, at 
the points rg defined by Scherk's surface. The 
value of 5<p (rg) vanishes — ► ±oo because 
Scherck's surface and the intermaterial dividing 
surface of our solution, defined by 6(j)(ri) = 0, co- 
incide in that limit. A convenient measure of the 
similarity of the two dividing surfaces, therefore, 
can be defined by computing 

h = J dr[</> A (r)-l/2] 2 <5(r-r s ) 

A measure /lsd for the LSD can be defined in 
the same manner. 
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Figure 5: Measures of the similarity between the 
intermaterial dividing surface defined by our cal- 
culation and that of Scherck's first surface (dia- 
monds, joined by solid line) and that defined by 
the linear stack of dislocations (squares, dotted 
line). Except near twist angles of ir/2, the latter 
is a better representation. 



A second means to compare the surfaces is to 
calculate the volume of the region which is en- 
closed between the two surfaces to be compared. 
This is easy to implement by a Monte Carlo in- 
tegration technique in which points of the unit 
cell are taken at random and checked to deter- 
mine whether the two surfaces agree, or not, in 
the assignment of the point to the A-rich re- 
gion. The relevant quantity is (5V)A/VN 1 / 2 a, 
the fraction of the volume SV for which there 
is disagreement, normalized by the area of the 
gram boundary V/A. The factor N 1 / 2 a ensures 
that this measure, which we denote ig, is dimen- 
sionless. 

These two measures are plotted in Fig. 5. 
They indicate that the LSD is a better represen- 
tation of the intermaterial dividing surface over 
almost the entire range of twist angles except, 
perhaps, quite close to ir/2. As the angle de- 
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Figure 6: Free energy per unit area, 7, of the 
twist grain boundary as a function of twist an- 
gle a. The circles show our results, the dashed 
line is a fit to them (extrapolated to a = 0) 
taking into account only data for a > 0.5. Also 
shown are the approximate free energies we have 
calculated from the expressions in Ref. 



121 for 



Scherck's first surface (solid line) and the linear 
array of dislocations (dotted line) which include 
only bending and compression energies. 



creases, both measures tend to the same limit, 
of course. (Recall the only difference between 
the two surfaces is a dilation of cos (a/2).) 

The grain boundary free energy as a function 
of the twist angle is shown in Fig. 6. The cir- 
cles show the results of our calculation. For val- 
ues of the twist angle greater than 0.5 radians, 
our results are accurate to within 5%, improving 
towards 1% as they approach it/ 2. The dashed 
line has been drawn through these values and ex- 
trapolated to zero. We have also included several 
other points for smaller values of a. The number 
of basis functions employed is insufficient for the 
grain boundary free energy to have converged to 
within 5%. Nevertheless, we have included them 
as they appear to indicate that the behavior of 



the free energy as the angle approaches zero may 
not be linear, the behavior expected if, at very 
low densities, dislocations repel one another fi"2|| . 
It has been argued, however, that at very low 
densities, dislocations attract one another [16]. 

One notes that the grain boundary free ener- 
gies are rather small: that of the grain boundary 
with twist angle of ir/2 is, at the same incompat- 
ibility, ^JV, somewhat less than half the energy 
of the boundary with tilt angle of 7r/2||. Per- 
haps this should not be too surprising. In the 
approximation noted earlier of treating the in- 
termaterial dividing surface as a surface of con- 
stant mean curvature, with an energy given by 



the Helfrich free energy[17], the grain boundary 
free energy would be identically zerofi~8|]. 

The free energy of a system with a twist grain 
boundary was calculated previously by Gido and 
Thomas p9|. They applied a version of the self 
consistent field to a brush of infinitely stretched 
chains anchored to a given saddle-shaped sur- 
face, and also carried out an independent calcu- 
lation based on work of Wang and Safran[17]. 
However they report their results in terms of 
the extensive free energy per chain in the region 
of Scherck's surface. This is not a uniquely de- 
fined quantity, nor is it the thermodynamic grain 
boundary free energy per unit area which we 
have calculated, so direct comparisons are pre- 
cluded. 

We have chosen to compare our results with 
those of Kamien and Lubensky[12|. We empha- 
size that the two calculations are rather different 
in principle. In the approach we have employed, 
the free energy of the block copolymer system is 
calculated directly, assuming nothing other than 
the applicability of self consistent field theory. In 
particular, we do not employ elasticity theory, or 
assume that displacements from a reference sys- 
tem without a grain boundary are small, etc. In 
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contrast, that of Ref. [12] assumes that the bulk 
system can be adequately described as a series 
of surfaces, and the energy of this system of sur- 
faces can be expanded assuming small displace- 
ments. A further difficulty which arises when 



ply be a mapping of the results of Ref. [12] to 



ours, and no independent comparison would be 
possible. Using these parameters, we have eval- 



applying the calculation of Ref. [12| to a block 
copolymer system is that the volumes on either 
side of the surfaces are undifferentiated, whereas 
in the block copolymer system, these volumes are 
filled with different monomers. Thus the sym- 
metry of the system considered by Kamien and 
Lubensky is not the same as that of ours. As a 
consequence, there are more elastic constants in 
an elastic description of a block copolymer lamel- 
lar phase than the two utilized by them in their 



description of liquid crystalline smectics[20[. 

Having acknowledged these caveats, we calcu- 
late the bending and compression contributions 
to the free energies of Scherck's first surface and 
of an LSD surface as given in Ref. [12]. As 
noted in the Introduction, a number of param- 
eters needed to evaluate these free energies are 
unknown, being inputs to the phenomenologi- 
cal theory. However we can provide some of 
these values from our work. Thus we calculate 
the lamellar spacing to be D = 1.5155 aN 1 / 2 , 
and the dimensionless compression modulus to 
be B = 3.01. The dimensionless bending modu- 
lus is unknown but can be estimated [^] to be 
K ~ 0.115. The only unknown which remains 
is the size of the "core region" , which provides a 
cutoff to the otherwise divergent integrals for the 
compression free energy estimated in |T^]. This 
can be estimated from the slope of the dotted 
line in Fig. 6. We obtain, thereby, a value of 
ps O.WlaN 1 / 2 , i.e., 25% of the IMDS spacing 
D/2. Of course, there is no reason that this core 
region should not depend on a, but were we to 
obtain the size of the core from our data at each 
value of the twist angle, the result would sim- 



uated the free energy given in Ref. [12] of the 
appropriate Scherck's surface and of the linear 
superposition of dislocations. These are shown 
in Fig. 6 as solid and dotted lines respectively. 
The LSD has a free energy closer to our result 
than does Scherck's surface, just as it is closer 
to our intermaterial dividing surface. Both ap- 
proximations underestimate the grain boundary 
free energy of the block copolymer system by a 
factor which increases with twist angle, and is 
about two at a = ir/2. 

4 Conclusions and outlook 

We have applied self consistent field theory to 
twist grain boundaries in block copolymer melts. 
Our calculation is more direct than earlier ones 
and provides greater information concerning the 
monomer densities throughout the volume. It 
also expresses the grain boundary free energy 
in terms of the directly measurable volume per 
chain and radius of gyration as opposed to elas- 
tic modulii of internal interfaces. The boundary 
free energy was obtained function of twist 
angle, and found to be quite small; smaller than 
kink grain boundaries of the same angle and in- 
compatibility. 

We have compared our results to previous phe- 
nomenological calculations to show that the in- 
termaterial dividing surface is not given by either 
surface, Scherck's first surface or the linear stack 
of dislocations, but that of the two, the latter is 
a better representation over most twist angles 
except near tt/2. 

We comment briefly on the other type of twist 
grain boundary which has been observed at small 
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twist angles || , the one consisting on a stack of 
lamellae which are twisted slightly and remain 
continuous. We have not investigated it is be- 
cause we failed to find an appropriate periodic 
boundary condition which does not contribute 
to the excess surface free energy. Although it is 
possible to calculate the contribution to the ex- 
cess surface free energy of any choice of boundary 
condition and then to subtract it from the total 
excess, leaving the desired grain boundary free 
energy, the procedure is tedious. However sim- 
ple examination of this boundary leads to the 
conclusion that the grain boundary free energy 
must be approximately twice that of a kink grain 
boundary. This is because the lamellae within 
the boundary and those far from it meet in what 
approximates a kink grain boundary. (Fig. 4 of 
Ref. |8| shows this nicely.) The free energy of a 
kink grain boundary grows as # 3 for small kink 
angle 9 ||] . Of course we do not know the rela- 
tion between the angle, 9, of this "effective" kink 
grain boundary and the twist angle a. Nonethe- 
less, if we assume that the relation is linear, then 
the twist grain boundary energy would grow as 
a 3 for small twist angles and would be favored 
over those we have modeled here, which would in 
fact be metastable but long-lived as their energy 
is small. This is in accord with the experimental 
result that both forms of boundary are observed 
at small twist angles. But at larger angles such 
grain boundaries would be disfavored compared 
to those considered in this paper. This is in ac- 
cord with the fact that they are not seen exper- 
imentally. 
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